Pred-O3, a web server to predict molecules, olfactory receptors and odor relationships.

The sense of smell is a biological process involving volatile molecules that interact with proteins called olfactory receptors to transmit a nervous message that allows the recognition of a perceived odor. However, the relationships between odorant molecules, olfactory receptors and odors (O3) are far from being well understood due to the combinatorial olfactory codes and large family of olfactory receptors. This is the reason why, based on 5802 odorant molecules and their annotations to 863 olfactory receptors (human) and 7029 odors and flavors annotations, a web server called Pred-O3 has been designed to provide insights into olfaction. Predictive models based on Artificial Intelligence have been developed allowing to suggest olfactory receptors and odors associated with a new molecule. In addition, based on the encoding of the odorant molecule's structure, physicochemical features related to odors and/or olfactory receptors are proposed. Finally, based on the structural models of the 98 olfactory receptors a systematic docking protocol can be applied and suggest if a molecule can bind or not to an olfactory receptor. Therefore, Pred-O3 is well suited to aid in the design of new odorant molecules and assist in fragrance research and sensory neuroscience. Pred-O3 is accessible at ' https://odor.rpbs.univ-paris-diderot.fr/'.

Exploring a Structural Data Mining Approach to Design Linkers for Head-to-Tail Peptide Cyclization.

Peptides have recently regained interest as therapeutic candidates, but their development remains confronted with several limitations including low bioavailability. Backbone head-to-tail cyclization, i.e., setting a covalent peptide bond linking the last amino acid with the first one, is one effective strategy of peptide-based drug design to stabilize the conformation of bioactive peptides while preserving peptide properties in terms of low toxicity, binding affinity, target selectivity, and preventing enzymatic degradation. Starting from an active peptide, it usually requires the design of a linker of a few amino acids to make it possible to cyclize the peptide, possibly preserving the conformation of the initial peptide and not affecting its activity. However, very little is known about the sequence-structure relationship requirements of designing linkers for peptide cyclization in a rational manner. Recently, we have shown that large-scale data-mining of available protein structures can lead to the precise identification of protein loop conformations, even from remote structural classes. Here, we transpose this approach to linkers, allowing head-to-tail peptide cyclization. First we show that given a linker sequence and the conformation of the linear peptide, it is possible to accurately predict the cyclized peptide conformation. Second, and more importantly, we show that it seems possible to elaborate on the information inferred from protein structures to propose effective candidate linker sequences constrained by length and amino acid composition, providing the first framework for the rational design of head-to-tail cyclization linkers. Finally, we illustrate this for two peptides using a limited set of amino-acids likely not to interfere with peptide function. For a linear peptide derived from Nrf2, the peptide cyclized starting from the experimental structure showed a 26-fold increase in the binding affinity. For urotensin II, a peptide already cyclized by a disulfide bond that exerts a broad array of biological activities, we were able, starting from models of the structure, to design a head-to-tail cyclized peptide, the first synthesized bicyclic 14-residue long urotensin II analogue, showing a retention of in vitro activity. Although preliminary, our results strongly suggest that such an approach has strong potential for cyclic peptide-based drug design.

PEP-FOLD4: a pH-dependent force field for peptide structure prediction in aqueous solution.

Accurate and fast structure prediction of peptides of less 40 amino acids in aqueous solution has many biological applications, but their conformations are pH- and salt concentration-dependent. In this work, we present PEP-FOLD4 which goes one step beyond many machine-learning approaches, such as AlphaFold2, TrRosetta and RaptorX. Adding the Debye-Hueckel formalism for charged-charged side chain interactions to a Mie formalism for all intramolecular (backbone and side chain) interactions, PEP-FOLD4, based on a coarse-grained representation of the peptides, performs as well as machine-learning methods on well-structured peptides, but displays significant improvements for poly-charged peptides. PEP-FOLD4 is available at http://bioserv.rpbs.univ-paris-diderot.fr/services/PEP-FOLD4. This server is free and there is no login requirement.

In Silico and In Vivo Studies of a Tumor-Penetrating and Interfering Peptide with Antitumoral Effect on Xenograft Models of Breast Cancer.

The combination of a tumor-penetrating peptide (TPP) with a peptide able to interfere with a given protein-protein interaction (IP) is a promising strategy with potential clinical application. Little is known about the impact of fusing a TPP with an IP, both in terms of internalization and functional effect. Here, we analyze these aspects in the context of breast cancer, targeting PP2A/SET interaction, using both in silico and in vivo approaches. Our results support the fact that state-of-the-art deep learning approaches developed for protein-peptide interaction modeling can reliably identify good candidate poses for the IP-TPP in interaction with the Neuropilin-1 receptor. The association of the IP with the TPP does not seem to affect the ability of the TPP to bind to Neuropilin-1. Molecular simulation results suggest that peptide IP-GG-LinTT1 in a cleaved form interacts with Neuropilin-1 in a more stable manner and has a more helical secondary structure than the cleaved IP-GG-iRGD. Surprisingly, in silico investigations also suggest that the non-cleaved TPPs can bind the Neuropilin-1 in a stable manner. The in vivo results using xenografts models show that both bifunctional peptides resulting from the combination of the IP and either LinTT1 or iRGD are effective against tumoral growth. The peptide iRGD-IP shows the highest stability to serum proteases degradation while having the same antitumoral effect as Lin TT1-IP, which is more sensitive to proteases degradation. Our results support the development of the TPP-IP strategy as therapeutic peptides against cancer.

Binding and Kinetic Analysis of Human Protein Phosphatase PP2A Interactions with Caspase 9 Protein and the Interfering Peptide C9h.

The serine/threonine phosphatase PP2A and the cysteine protease Caspase 9 are two proteins involved in physiological and pathological processes, including cancer and apoptosis. We previously demonstrated the interaction between Caspase 9 and PP2A and identified the C9h peptide, corresponding to the binding site of Caspase 9 to PP2A. This interfering peptide can modulate Caspase 9/PP2A interaction leading to a strong therapeutic effect in vitro and in vivo in mouse models of tumor progression. In this manuscript, we investigate (I) the peptide binding to PP2A combining docking with molecular dynamics and (II) the secondary structure of the peptide using CD spectroscopy. Additionally, we compare the binding affinity, using biolayer interferometry, of the wild-type protein PP2A with Caspase 9 and vice versa to that observed between the PP2A protein and the interfering peptide C9h. This result strongly encourages the use of peptides as new therapeutics against cancer, as shown for the C9h peptide already in clinical trial.

SeamDock: An Interactive and Collaborative Online Docking Resource to Assist Small Compound Molecular Docking.

In silico assessment of protein receptor interactions with small ligands is now part of the standard pipeline for drug discovery, and numerous tools and protocols have been developed for this purpose. With the SeamDock web server, we propose a new approach to facilitate access to small molecule docking for nonspecialists, including students. The SeamDock online service integrates different docking tools in a common framework that allows ligand global and/or local docking and a hierarchical approach combining the two for easy interaction site identification. This service does not require advanced computer knowledge, and it works without the installation of any programs with the exception of a common web browser. The use of the Seamless framework linking the RPBS calculation server to the user's browser allows the user to navigate smoothly and interactively on the SeamDock web page. A major effort has been put into the 3D visualization of ligand, receptor, and docking poses and their interactions with the receptor. The advanced visualization features combined with the seamless library allow a user to share with an unlimited number of collaborators, a docking session, and its full visualization states. As a result, SeamDock can be seen as a free, simple, didactic, evolving online docking resource best suited for education and training.

Molecular dynamics simulations reveal statistics and microscopic mechanisms of water permeation in membrane-embedded artificial water channel nanoconstructs.

Understanding water transport mechanisms at the nanoscale level remains a challenge for theoretical chemical physics. Major advances in chemical synthesis have allowed us to discover new artificial water channels, rivaling with or even surpassing water conductance and selectivity of natural protein channels. In order to interpret experimental features and understand microscopic determinants for performance improvements, numerical approaches based on all-atom molecular dynamics simulations and enhanced sampling methods have been proposed. In this study, we quantify the influence of microscopic observables, such as channel radius and hydrogen bond connectivity, and of meso-scale features, such as the size of self-assembly blocks, on the permeation rate of a self-assembled nanocrystal-like artificial water channel. Although the absolute permeation rate extrapolated from these simulations is overestimated by one order of magnitude compared to the experimental measurement, the detailed analysis of several observed conductive patterns in large assemblies opens new pathways to scalable membranes with enhanced water conductance for the future design.

InterEvDock3: a combined template-based and free docking server with increased performance through explicit modeling of complex homologs and integration of covariation-based contact maps.

The InterEvDock3 protein docking server exploits the constraints of evolution by multiple means to generate structural models of protein assemblies. The server takes as input either several sequences or 3D structures of proteins known to interact. It returns a set of 10 consensus candidate complexes, together with interface predictions to guide further experimental validation interactively. Three key novelties were implemented in InterEvDock3 to help obtain more reliable models: users can (i) generate template-based structural models of assemblies using close and remote homologs of known 3D structure, detected through an automated search protocol, (ii) select the assembly models most consistent with contact maps from external methods that implement covariation-based contact prediction with or without deep learning and (iii) exploit a novel coevolution-based scoring scheme at atomic level, which leads to significantly higher free docking success rates. The performance of the server was validated on two large free docking benchmark databases, containing respectively 230 unbound targets (Weng dataset) and 812 models of unbound targets (PPI4DOCK dataset). Its effectiveness has also been proven on a number of challenging examples. The InterEvDock3 web interface is available at http://bioserv.rpbs.univ-paris-diderot.fr/services/InterEvDock3/.

Coarse-grained and atomic resolution biomolecular docking with the ATTRACT approach.

The ATTRACT protein-protein docking program has been employed to predict protein-protein complex structures in CAPRI rounds 38-45. For 11 out of 16 targets acceptable or better quality solutions have been submitted (~70%). It includes also several cases of peptide-protein docking and the successful prediction of the geometry of carbohydrate-protein interactions. The option of combining rigid body minimization and simultaneous optimization in collective degrees of freedom based on elastic network modes was employed and systematically evaluated. Application to a large benchmark set indicates a modest improvement in docking performance compared to rigid docking. Possible further improvements of the docking approach in particular at the scoring and the flexible refinement steps are discussed.

Glutathionylation primes soluble glyceraldehyde-3-phosphate dehydrogenase for late collapse into insoluble aggregates.

Protein aggregation is a complex physiological process, primarily determined by stress-related factors revealing the hidden aggregation propensity of proteins that otherwise are fully soluble. Here we report a mechanism by which glycolytic glyceraldehyde-3-phosphate dehydrogenase of Arabidopsis thaliana (AtGAPC1) is primed to form insoluble aggregates by the glutathionylation of its catalytic cysteine (Cys149). Following a lag phase, glutathionylated AtGAPC1 initiates a self-aggregation process resulting in the formation of branched chains of globular particles made of partially misfolded and totally inactive proteins. GSH molecules within AtGAPC1 active sites are suggested to provide the initial destabilizing signal. The following removal of glutathione by the formation of an intramolecular disulfide bond between Cys149 and Cys153 reinforces the aggregation process. Physiological reductases, thioredoxins and glutaredoxins, could not dissolve AtGAPC1 aggregates but could efficiently contrast their growth. Besides acting as a protective mechanism against overoxidation, S-glutathionylation of AtGAPC1 triggers an unexpected aggregation pathway with completely different and still unexplored physiological implications.

DaReUS-Loop: a web server to model multiple loops in homology models.

Loop regions in protein structures often have crucial roles, and they are much more variable in sequence and structure than other regions. In homology modeling, this leads to larger deviations from the homologous templates, and loop modeling of homology models remains an open problem. To address this issue, we have previously developed the DaReUS-Loop protocol, leading to significant improvement over existing methods. Here, a DaReUS-Loop web server is presented, providing an automated platform for modeling or remodeling loops in the context of homology models. This is the first web server accepting a protein with up to 20 loop regions, and modeling them all in parallel. It also provides a prediction confidence level that corresponds to the expected accuracy of the loops. DaReUS-Loop facilitates the analysis of the results through its interactive graphical interface and is freely available at http://bioserv.rpbs.univ-paris-diderot.fr/services/DaReUS-Loop/.

Water permeation across artificial I-quartet membrane channels: from structure to disorder.

Artificial water channels (AWCs) have been designed for water transport across membranes with the aim to mimic the high water permeability observed for biological systems such as aquaporins (∼108-109 water molecules per s per channel), as well as their selectivity to reject ion permeation at the same time. Recent works on designed self-assembling alkylureido-ethylimidazole compounds forming imidazole-quartet channels (I-quartets), have shown both high water permeability and total ionic-rejection. I-quartets are thus promising candidates for further development of AWCs. However, the molecular mechanism of water permeation as well as I-quartet organization and stability in a membrane environment need to be fully understood to guide their optimal design. Here, we use a wide range of all-atom molecular dynamics (MD) simulations and their analysis to understand the structure/activity relationships of the I-quartet channels. Four different types with varying alkyl chain length or chirality have been studied in a complex fully hydrated lipid bilayer environment at both microsecond and nanosecond scale. Microsecond simulations show two distinct behaviors; (i) two out of four systems maintain chiral dipolar oriented water wires, but also undergo a strong reorganization of the crystal shape, (ii) the two other I-quartet channels completely lose the initial organization, nonetheless keeping a water transport activity. Short MD simulations with higher time resolution were conducted to characterize the dynamic properties of water molecules in these model channels and provided a detailed hypothesis on the molecular mechanism of water permeation. The ordered confined water was characterized with quantitative measures of hydrogen-bond life-time and single particle dynamics, showing variability among I-quartet channels. We will further discuss the underlying assumptions, currently based on self-aggregation simulations and crystal patches embedded in lipid bilayer simulations and attempt to describe possible alternative approaches to computationally capture the water permeation mechanism and the self-assembly process of these AWCs.

Oriented chiral water wires in artificial transmembrane channels.

Aquaporins (AQPs) feature highly selective water transport through cell membranes, where the dipolar orientation of structured water wires spanning the AQP pore is of considerable importance for the selective translocation of water over ions. We recently discovered that water permeability through artificial water channels formed by stacked imidazole I-quartet superstructures increases when the channel water molecules are highly organized. Correlating water structure with molecular transport is essential for understanding the underlying mechanisms of (fast) water translocation and channel selectivity. Chirality adds another factor enabling unique dipolar oriented water structures. We show that water molecules exhibit a dipolar oriented wire structure within chiral I-quartet water channels both in the solid state and embedded in supported lipid bilayer membranes (SLBs). X-ray single-crystal structures show that crystallographic water wires exhibit dipolar orientation, which is unique for chiral I-quartets. The integration of I-quartets into SLBs was monitored with a quartz crystal microbalance with dissipation, quantizing the amount of channel water molecules. Nonlinear sum-frequency generation vibrational spectroscopy demonstrates the first experimental observation of dipolar oriented water structures within artificial water channels inserted in bilayer membranes. Confirmation of the ordered confined water is obtained via molecular simulations, which provide quantitative measures of hydrogen bond strength, connectivity, and the stability of their dipolar alignment in a membrane environment. Together, uncovering the interplay between the dipolar aligned water structure and water transport through the self-assembled I-quartets is critical to understanding the behavior of natural membrane channels and will accelerate the systematic discovery for developing artificial water channels for water desalting.

Simulation of Ligand Binding to Membrane Proteins.

Membrane proteins are involved in a large variety of functions. Most of these protein functions are regulated by ligand binding with diverse modes of action: agonists, partial agonists, antagonists, and allosteric modulators, potentiators and inhibitors. From the pharmacological point of view, membrane proteins are one if not the major target for drug development. However, experimental structure determination of membrane proteins in complex or in free form still represents a great challenge. Molecular dynamics (MD) simulations commonly reach the microsecond scale on membrane systems. This numerical tool is mature enough to predict and add molecular details on the different ligand-binding modes. In the present chapter, I will present the different steps to design, simulate, and analyze a MD simulation system containing a protein embedded in a membrane and surrounded by water and ligand. As an illustration, the simulation of the ligand-gated ion channel γ-aminobutyric acid type A receptor (GABAAR) surrounded by one of its allosteric potentiators, bromoform, will be presented and discussed.

String method solution of the gating pathways for a pentameric ligand-gated ion channel.

Pentameric ligand-gated ion channels control synaptic neurotransmission by converting chemical signals into electrical signals. Agonist binding leads to rapid signal transduction via an allosteric mechanism, where global protein conformational changes open a pore across the nerve cell membrane. We use all-atom molecular dynamics with a swarm-based string method to solve for the minimum free-energy gating pathways of the proton-activated bacterial GLIC channel. We describe stable wetted/open and dewetted/closed states, and uncover conformational changes in the agonist-binding extracellular domain, ion-conducting transmembrane domain, and gating interface that control communication between these domains. Transition analysis is used to compute free-energy surfaces that suggest allosteric pathways; stabilization with pH; and intermediates, including states that facilitate channel closing in the presence of an agonist. We describe a switching mechanism that senses proton binding by marked reorganization of subunit interface, altering the packing of ß-sheets to induce changes that lead to asynchronous pore-lining M2 helix movements. These results provide molecular details of GLIC gating and insight into the allosteric mechanisms for the superfamily of pentameric ligand-gated channels.

Salt-Excluding Artificial Water Channels Exhibiting Enhanced Dipolar Water and Proton Translocation.

Aquaporins (AQPs) are biological water channels known for fast water transport (∼10(8)-10(9) molecules/s/channel) with ion exclusion. Few synthetic channels have been designed to mimic this high water permeability, and none reject ions at a significant level. Selective water translocation has previously been shown to depend on water-wires spanning the AQP pore that reverse their orientation, combined with correlated channel motions. No quantitative correlation between the dipolar orientation of the water-wires and their effects on water and proton translocation has been reported. Here, we use complementary X-ray structural data, bilayer transport experiments, and molecular dynamics (MD) simulations to gain key insights and quantify transport. We report artificial imidazole-quartet water channels with 2.6 Špores, similar to AQP channels, that encapsulate oriented dipolar water-wires in a confined chiral conduit. These channels are able to transport ∼10(6) water molecules/s, which is within 2 orders of magnitude of AQPs' rates, and reject all ions except protons. The proton conductance is high (∼5 H(+)/s/channel) and approximately half that of the M2 proton channel at neutral pH. Chirality is a key feature influencing channel efficiency.

Sites of Anesthetic Inhibitory Action on a Cationic Ligand-Gated Ion Channel.

Pentameric ligand-gated ion channels have been identified as the principal target of general anesthetics (GA), whose molecular mechanism of action remains poorly understood. Bacterial homologs, such as the Gloeobacter violaceus receptor (GLIC), have been shown to be valid functional models of GA action. The GA bromoform inhibits GLIC at submillimolar concentration. We characterize bromoform binding by crystallography and molecular dynamics (MD) simulations. GLIC's open form structure identified three intra-subunit binding sites. We crystallized the locally closed form with an additional bromoform molecule in the channel pore. We systematically compare binding with the multiple potential sites of allosteric channel regulation in the open, locally closed, and resting forms. MD simulations reveal differential exchange pathways between sites from one form to the other. GAs predominantly access the receptor from the lipid bilayer in all cases. Differential binding affinity among the channel forms is observed; the pore site markedly stabilizes the inactive versus active state.